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1.  INTRODUCTION 


The  ability  to  sense  a  penetration  event  and  determine  the  point  of  impact  on  a  two- 
dimensional  surface  is  a  necessary  prerequisite  for  many  existing  and  future  armor  systems. 
The  degree  of  sophistication  of  these  measurements  covers  a  wide  range  and  is  mandated  by 
the  testing  or  systems  requirements.  For  some  studies,  such  as  evaluation  of  aiming  and 
projectile  accuracy,  a  ”post-evenr  determination  of  impact  point  is  suffident  This  type  of 
analysis  can  be  as  mdimentary  as  finding  where  the  holes  are  in  a  target  The  time  of  the 
penetration  event  needs  to  be  accurately  and  rapidly  detennined  if  the  dynamics  of  the 
projectile-armor  interaction  are  to  be  investigated  by  flash  x-ray.  high-speed  optical  cameras, 
or  other  analysis  equipment  For  high  resolution  analysis,  the  location  of  the  penetration 
event  may  also  need  to  be  determined. 

The  break-wire  saeen  has  been  the  standard  tool  for  sensing  penetration  events  and 
measuring  their  location  in  real-time.  A  typicsti  embodiment  of  this  apparatus  indudes  a  series 
of  parallel,  isolated  electrical  conductors  arranged  in  a  plane  which  is  normal  to  the  path  of  the 
projectile.  For  a  two-dimensional  determination  of  location,  a  second  array  is  induded  where 
the  conductors  are  perpendicular  with  respect  to  the  conductors  in  the  first  array.  As  a 
projectile  passes  through  this  grid,  it  breaks  soma  of  the  conductors,  causing  a  loss  of 
electrical  contact  The  time  of  the  initial  penetration  event  can  be  denoted  by  feeding  all  the 
signals  from  the  individual  conductors  into  an  AND  gate  and  sensing  the  change  of  state  of 
the  output  A  subsequent  poll  of  the  states  of  the  individual  conductors  reveals  the  location  of 
the  penetration  event 

Advantages  of  the  break-wire  method  include  the  overall  simplicity  of  the  concept,  its 
immunity  to  noise,  and  its  intrinsically  digital  nature.  However,  in  spite  of  the  simplicity  of  the 
concept  the  standard  application  of  this  technique  can  present  substantial  engineering 
difficulties.  The  requirement  of  numerous  isolated  conductors  and  separate  circuitry  for  each 
wire,  or  multiplexers,  to  monitor  the  state  of  the  screen  can  make  the  task  of  installing  and 
maintaining  such  a  system  tedious  or  even  impractical.  To  ensure  that  a  projectile  will  break 
at  least  one  conductor,  the  spacing  between  adjacent  conductors  can  not  exceed  the 
projectile  diameter.  Therefore,  complexity  grows  rapidly  when  the  impact  site  of  small  objects 
must  be  established  over  a  large  area 
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The  present  study  seeks  development  of  an  alternative  method  for  measuring  Impact  sites 
on  a  two-dimensional  surface.  A  proof-of-prindpie  device  resulted  which  can  easily  be  scaled 
up  to  monitor  large  surfaces  and  has  the  ability  to  determine  positions  before  the  impacting 
ballistic  projectile  can  travel  a  significant  distance.  The  solution  offered  here  rests  upon 
creating  a  charge  flow  from  the  impact  point  through  an  electrically  resistive  medium  to  the 
extremities  of  the  sensing  surface.  A  measurement  of  the  relative  rate  of  flow  to  each  sensor 
edge  indicates  the  impact  location.  In  principle,  this  type  of  sensing  surface  can  be  expanded 
to  any  reasonable  size  without  changing  the  nature  of  its  operation.  The  determination  of  the 
impact  location  in  the  proof-of-principle  device  takes  a  few  tens  of  microseconds,  a  figure 
which  can  doubtiess  be  reduced  substantially.  Rnaiiy,  the  device  is  inexpensive,  simple,  and 
can  be  fabricated  from  readily  available  materials  and  components. 

2.  EXPERIMENTAL  PROCEDURE 

2.1  System  Hardware. 

2.1.1  Position  Sensing  Transducer.  The  active  surface  of  the  position  sensor  consists  of 
an  electrically  resistive  square  sheet  with  four  charge  collection  tabs  fixed  along  the  outer 
edges.  The  geometry  of  the  active  surface  and  collection  tabs  is  illustrated  in  Figure  1 .  When 
electrical  charge  arrives  on  this  surface,  it  raises  the  local  potentir.!,  causing  the  charge  to  flow 
outward  to  the  collection  tabs.  The  resistance  of  the  sheet  is  set  so  that  the  rate  at  which 
charge  is  collected  by  the  four  tabs  depends  on  their  distance  to  the  point  where  the  charge  is 
initially  deposited.  The  closer  a  tab  is  to  the  point  of  charge  deposition,  the  lower  the  effective 
path  resistance  to  the  tab  and  the  higher  the  rate  of  charge  transfer  and  collection.  By 
com^jaring  the  relative  amounts  of  charge  collected  by  the  tabs  on  opposing  sides  of  the 
sensing  surface  (e.g.,  X,  vs.  Xj  and  vs.  Yj),  the  position  of  the  point  of  charge  deposition 
along  each  axis  can  be  deduced. 

The  resistive  sheet  material  is  prepared  by  impregnating  vellum  paper  with  fine  grain 
graphite  particles.  This  recipe  arose  from  tests  of  several  fabrication  schemes,  mostly 
involving  various  means  of  spraying  and  soaking  of  graphite  in  a  volatile  solvent  The  simple 
rubbing  of  the  graphite  powder  into  the  paper  surface  is  the  easiest  and  yields  the  most 
uniform  sheet  resistivities.  Although  beyond  the  scope  of  this  report,  it  may  be  possible  to 
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optimize  the  perfotmance  of  the  position  sensor  by  adjusting  the  sheet  resistivity  of  the 
impregnated  particulate  material  by  the  addition  of  conducting/insulating  materials  to  the 
graphite  base  and  by  adjusting  the  distribution  of  resistivity  on  the  active  surface.  All  of  the 
experimental  results  presented  in  this  report  were  obtained  using  a  20-cm  by  20-cm  sensing 
surface.  The  collection  tabs  consisted  of  1 0-cm-long  blocks  which  were  covered  with  copper- 
foil  tape  and  mechanically  conriS's.ted  to  the  edges  of  the  resistive  sheet  to  yield  an 
insignificant  electrical  contact  resistar^e.  The  nominal  resistance  on  the  sensing  surface  was 
150  ohms  between  opposing  tabs  and  125  ohms  between  normal  tabs. 

There  are  a  variety  of  potential  methods  for  depositing  charge  to  the  position  sensor’s 
active  surface.  For  field  applications,  a  small  capacitor  can  be  connected  to  a  thin  metallic 
sheet  which  closely  covers,  but  does  not  actually  touch,  the  active  surface.  When  impacted 
by  an  object,  the  metallic  sheet  deforms  and  deposits  the  charge  from  the  capacitor  at  the 
point  of  impact  With  an  appropriate  increase  in  surface  resistivity,  piezoelectric  plastic  films 
might  be  utilized  in  a  similar  manner  with  the  added  advantage  of  making  the  discharge 
capacitor  unnecessary. 
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For  laboratory  testing  and  calibration,  it  is  important  to  repeatedly  deposit  charge  on  the 
active  surface  at  known  locations.  To  facilitate  such  testing,  a  metallic  cover  plate  which  had 
an  array  of  small  holes  spanning  the  entire  active  surface  of  the  position  detector  was  used. 

By  passing  a  fine  metallic  probe  through  these  holes,  the  cover  plate  could  be  shorted  to  the 
active  surface  and  the  charge  deposited  at  a  well  defined  location.  The  charge  was  stored  in 
a  220-pF  capacitor  charged  to  -15  V.  A  5,100-ohm  resistor  was  connected  between  the 
capacitor  and  cover  plate  to  regulate  the  discharge  oirrent  The  charge  collection  tabs  were 
coupled  through  0.1 -^F  capacitors  to  the  signal  processing  electronics.  This  AC  coupling 
offered  a  high  resistance  to  60-cycie  noise.  A  large  resistor  was  also  connected  between  the 
active  surface  and  ground  to  prevent  the  undesirable  buildup  of  deposited  charge. 

2.A  .2  Signal  Processing.  A  significant  feature  of  the  position  sensor’s  signal  processing 
unit  is  the  relative  insensitivity  to  the  actual  amount  of  charge  which  is  deposited  to  the  active 
surface.  This  broad  range  of  acceptable  charge  deposition  is  achieved  by  measuring  the 
rates  at  which  charge  is  collected  rather  than  the  total  amount  of  charge.  An  added  benefit  of 
this  approach  is  the  enhanced  speed  made  possible  by  a  measurement  strategy  which 
operates  while  the  charge  is  being  transferred  rather  than  having  to  wait  until  all  the  charge 
has  been  collected.  An  overview  of  the  signal  processing  unit  is  illustrated  in  the  block 
diagram  of  Figure  2.  Detailed  schematics  and  descriptions  of  all  electronic  circuitry  are  in 
Appendix  A. 

The  first  step  in  measuring  the  rate  of  charge  transfer  is  to  convert  the  charge  collected  at 
the  four  tabs  into  amplified  voltage  signals.  This  task  is  performed  by  the  charge  amplifiers. 
When  the  sum  of  the  voltage  signals  for  each  pair  of  opposing  tabs  reaches  a  predetermined 
threshold  level,  the  individual  voltage  signals  and  their  sum  are  captured  by  sample  and  hold 
circuits.  An  analog  division  circuit  then  determines  the  value  of  the  quotient  fab  2  signal/(tab  1 
signal-htab  2  signal).  This  representation  for  combining  the  signals  from  opposing  tabs 
provides  results  that  are  a  linear  function  of  the  charge  collected  at  each  tab  between  the  well 
defined  limits  of  0  and  1 .  The  seemingly  simpler  representation  tab  1  signaL'tab  2  signal 
yields  a  distinctly  nonlinear  result  between  the  ill-defined  limits  of  0  and  <»>.  A  voltage  amplifier 
which  can  vary  the  DC  offset  as  well  as  the  amplitude  of  the  quotient  signal  is  used.  This 
allows  the  output  signal  to  be  set  at  a  convenient  level  for  viewing  by  cscilloscope  or 
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Figure  2.  Block  Diagram  of  Position  Sensor  Signal  Processing  Unit. 

other  analog  display  instrumentation.  In  addition,  the  voltage  amplifier  expands  the  range  of 
possible  analog  output  signals  to  completely  cover  the  measurement  range  of  subsequent 
analog  to  digital  converters  in  order  to  achieve  maximum  conversion  resolution.  These 
digitized  signals  are  readily  accessible  for  computer-controlled  calibration  and  display. 

2.2  System  Software. 

2.2.1  Software  Overview.  In  order  for  this  position  sensor  to  be  integrated  into  a  decision¬ 
making  or  control  environment,  its  output  signals  should  be  converted  to  digital  form.  Analog 
to  digital  (A/D)  conversion  of  the  voltage  amplifier  output  signals  (Figure  2)  is  performed  with 
the  assistance  of  computer  control.  The  use  of  computer  control  requires  associated  software 
drivers.  For  system  software  deveiopment.  an  Apple  lie  (Apple  Computer  Inc..  Cupertino.  CA) 
computer  is  used  as  the  controller  with  high  level  programs  written  in  Applesoft  BASIC  and 
machine  language  routines  written  for  the  6502  microprocessor.  The  1-MHz  clock  frequency 
and  interpreted  high  level  language  do  not  yield  a  particularly  fast  system  by  today’s 
standards.  However,  this  setup  is  more  than  adequate  for  determining  system  feasibility  and 
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developing  software  logic.  Digital  conversions  are  performed  by  an  8-bit  ADC0816  (National 
Semiconductor  Corporation,  Santa  Clara,  CA)  monolithic  analog-to-digital  converter  controlled 
by  the  microprocessor. 

The  software  that  controls  the  data  acquisition  system  can  be  divided  into  the  overlapping 
categories  of  system  calibration,  data  acquisition,  data  correction,  and  data  display.  A 
schematic  illustrating  the  interaction  between  the  software  units  that  perfonn  these  tasks  is 
shown  in  Rgure  3.  High  level  language  code  is  denoted  by  boxes  with  rounded  comers  while 
data  files  are  shown  as  rectangles.  Machine  code  that  is  accessed  by  a  high  level  language 
routine  is  indicated  by  a  dashed  box. 

2.2.2  Calibration  Software.  An  added  advantage  of  computer  control  is  the  potential  for 
sophisticated  system  calibration  and  data  correction.  The  ideal  output  characteristics  for  one 
dimension  of  a  two-dimensional  position  sensor  are  illustrated  in  Rgure  4.  In  this  example, 
the  determined  X  position  is  plotted  as  a  function  of  the  actuai  X  and  Y  positions,  X,  and  Y,. 
For  the  ideal  case,  the  determined  X  position  is  linearly  proportional  to  X,  and  is  not  a  function 
of  Y,.  The  other  dimension  would  have  similar  characteristics.  That  is,  the  determined  Y 
position  would  be  independent  of  X.  and  linearly  proportionaJ  to  Y,. 

Figure  5  illustrates  the  actual  raw  output  characteristics  for  the  position  sensor  in  the  X 
direction.  In  order  to  realize  maximum  resolution  and  accuracy,  the  zero  offset  and  scale  of 
the  voltage  amplifiers  is  first  adjusted  so  that  the  raw  X  values,  X,,  cover  roughly  the  entire 
range  of  allowe'J  values.  The  data  are  then  obtained  by  measuring  X,  values  for  each 
element  in  the  grid  of  known  actuai  positions  (X,,YJ  across  the  active  surface  of  the  position 
sensor.  In  contrast  to  the  ideal  characteristics  shown  in  Figure  4,  the  actual  output 
characteristics  are  neither  independent  of  Y,  nor  particularly  linear  with  respect  to  X,. 
Therefore,  any  data  correction  scheme  must  treat  a  corrected  value  for  the  X  position,  X^,  to 
be  a  function  of  both  the  raw  X  and  Y  data  values: 

Xc  =  f{W  (1) 
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Figure  4.  Ideal  Position  Sensor  X  Axis  Output  Characteristics. 
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Figure  5.  Uncorrect  Position  Sensor  X  Axis  Output  Characteristics. 

Similarty,  for  the  corrected  value  of  the  Y  position,  Y^: 

Yo  =  g(X,Y,).  (2) 

At  minimum  and  maximum  values  for  Y,.  the  raw  X,  values  for  centrally  located  actual 
positions,  X.,  are  pinned  to  a  median  value  by  the  presence  of  the  two  charge  collection  tabs 
which  are  simultaneously  measuring  the  Y  position.  Thus,  instead  of  defining  an  inclined 
ramp  in  position  space,  the  position  sensor's  raw  output  characteristics  define  a  curved 
saddle-like  shape  which  clearly  illustrates  the  need  for  data  correction  processing. 

The  plots  in  Figures  4  and  5  consider  the  output  characteristics  of  only  a  single  dimension 
of  the  two-dimensional  position  sensor.  In  order  to  simultaneously  represent  the  transducer 
response  of  both  axes,  the  plots  of  Figures  6  and  7  are  convenient.  Rgure  6  illustrates  the 
output  characteristics  from  an  idealized  two-dimensional  position  sensor.  The  determined 
positions  from  a  regular  array  of  actual  positions  {X,,YJ  generate  a  linear  grid.  By  contrast 
Figure  7  shows  that  only  the  centrally  located  positions  of  the  position  sensor  yield  an 
approximately  regular  grid  of  raw  data  values.  A  major  source  of  this  distortion  is  the 
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presence  of  the  charge  collection  tabs.  When  charge  is  deposited  in  a  comer  position,  the 
charge  collection  tabs  on  the  adjacent  edges  offer  a  low  resistance  path  to  the  other  side  of 
the  active  surface.  The  apparent  point  of  deposition  is  therefore  located  at  a  more  central 
position.  One  obvious  remedy  for  this  problem  would  be  to  reduce  the  length  of  the  charge 
collection  tabs.  However,  the  use  of  overly  short  tabs  can  lead  to  an  increased  number  of 
positions  which  yield  nonunique  raw  data  values.  Future  efforts  devoted  to  refining  this 
position  detector  should  include  studies  to  determine  the  optimal  charge  collection  tab 
geometries.  At  any  rats,  once  again  the  need  for  data  conection  processing  is  clear. 

A  first  step  in  calibrating  the  position  sensor  is  to  obtain  raw  caiibration  data  by  depositing 
charge  at  a  matrix  of  known  locations  on  the  sensor  surface.  This  calibration  data  is  of  the 
type  shown  in  Figure  5.  Each  value  of  raw  position  output  from  the  sensor  is  normalized  and 
rounded  off  to  integral  values  from  0  to  39.  In  this  representation,  the  smallest  increment  of 
position  corresponds  to  1/40th  of  full  scale.  This  scaling  is  performed  to  match  the  low 
resolution  graphics  capability  of  the  Apple  He.  whidi  will  eventually  also  serve  as  the  display 
device. 

To  measure  the  response  of  the  position  detector  at  all  1 ,600  locations  on  a  40  by  40  grid 
would  be  excessively  tedious.  Therefore,  data  are  acquired  on  a  more  manageable  17  by 
1 7  grid  which  is  stored  in  a  data  file.  Maximum  resolution  for  the  corrected  position  value 
ultimately  to  come  from  this  device  (Xc,YJ  still  requires  one  to  know  the  raw  data  (X,.Y,)  on  a 
full  40  by  40  grid  of  points.  Thus,  the  data  file  for  the  17  by  17  array  of  experimentally 
measured  raw  data  is  run  through  another  high  level  routine  which  interpolates  the  data  to  a 
40  by  40  array  of  points  and  stores  this  information  in  an  additional  file.  The  errors  introduced 
by  the  interpolation  procedure  are  insignificant 

The  array  is  then  inverted,  and  a  look-up  table  of  corrected  positions  associated  with  each 
possible  raw  data  pair  is  stored  in  a  calibration  data  file — a  file  which  is  indexed  by  its  raw 
data  values.  The  look-up  table,  which  in  operation  is  stored  in  computer  memory,  is  intended 
to  form  a  bridge  between  the  raw  and  corrected  data  sets.  When  the  device  responds  to  an 
impact  by  recording  a  raw  data  position  (X,.Y,),  a  high  level  routine  scans  the  memory  for  the 
address  indexed  by  this  pair  and  reports  the  corresponding  corrected  position  (X^.YJ  stored 
therein. 
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The  goal  of  this  effort  and  indeed  the  measure  of  this  device’s  worth  is  that  the  corrected 
position  shouid  be  equal  to  the  actual  position,  (X.,YJ.  The  high  ievel  routines  which 
interpoiate  the  calibration  raw  data  matrix  and  generate  the  look-up  table  are  listed  in 
Appendices  B-1  and  B-2. 

22.3  Data  Acquisition.  Correction,  and  Oisii^ay  Routine.  An  additional  high  level  routine, 
one  dedicated  to  controlling  actual  impact  data  acquisition,  is  listed  in  Appendix  B-3  along  with 
a  detailed  accounting  of  the  associated  memory  configuration.  This  routine  acquires  raw  data, 
uses  the  calibration  look-up  table  to  correct  the  data,  and  then  displays  the  result  To 
increase  the  speed  of  the  data  acquisition  process,  this  high  level  routine  calls  a  machine 
language  routine  to  control  the  A/0  chip  for  the  actual  data  acquisition.  The  assembler  listing 
of  this  machine  code  routine  is  also  listed  in  Appendix  B-4. 

After  the  machine  language  A/D  driver  is  called  and  the  raw  data  is  scaled  and  corrected, 
the  corrected  position  is  graphically  indicated  on  the  computer  monitor  by  the  intersection  of 
two  cross  hairs.  With  the  output  signal  from  the  position  sensor  available  in  both  analog  and 
digital  form,  a  wide  range  of  display,  control,  and  decision-making  procedures  could  now  be 
implemented. 

2.3  System  Test  and  Evaluation. 

2.3.1  Hardware  Testing.  When  charge  is  deposited  to  the  resistive  active  surface  of  the 
position  sensor,  it  immediately  begins  to  flow  to  ttie  array  of  tabs  along  the  edges.  The  rate  at 
which  charge  is  collected  at  each  tab  depends  on  the  electrical  resistance,  or  distance, 
between  the  deposition  point  and  the  tab  in  question.  Once  collected,  the  charge  is 
processed  by  a  charge  amplifier  (Rgure  A-1),  the  output  wavefonn  of  which  is  dependent  on 
the  position  of  charge  deposition.  Figure  8  shows  the  output  from  two  charge  amplifiers 
connected  to  opposing  tabs  along  one  dimension  of  the  position  sensor  for  a  series  of  charge 
deposition  locations.  An  illustration  of  the  relative  positions  of  the  tabs  and  charge  deposition 
point  is  shown  on  the  left  with  the  associated  oscillogram  of  the  charge  amplifier  outputs  on 
the  right  For  clarity,  the  baselines  of  the  two  oscilloscope  traces  have  been  offset  The 
scales  for  all  the  osdilograms  are  as  noted  in  the  upper  figure.  The  oscilloscope  is 
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pretriggered  to  allow  the  baseline  levels  to  be  recorded,  in  each  case,  the  actual  time  of 
triggering  is  at  or  near  the  first  indication  of  activity  above  the  base  level. 

When  charge  is  deposited  on  the  left  portion  of  the  sheet  the  distance  to  the  ^  tab  is 
much  less  than  the  distance  to  the  S  tab.  causing  the  effective  resistance  to  the  i4  tab  to  also 
be  less.  Thus,  the  rate  of  charge  collection  at  toe  ^  tab  should  be  greater  than  at  toe  8  tab, 
an  inference  corroborated  by  toe  greater  slope  of  trace  A  in  the  associated  oscillogram.  As 
the  charge  deposition  point  is  moved  to  the  right  the  rata  of  charge  collection  at  tab  A 
decreases  while  the  charge  collection  rate  at  tab  S  increases.  When  charge  is  deposited  in 
toe  middle  of  the  sheet  toe  collection  rate  at  the  two  tabs  is  approximately  equal  as  indicated 
by  toe  equal  slopes  of  the  two  charge  amplifier  output  signals.  As  toe  point  of  charge 
deposition  is  moved  farther  to  the  right  the  trend  continues  with  greater  charge  collection 
rates  and  amplifier  signal  slopes  being  noted  for  tab  S. 

The  output  from  the  charge  amplifiers  is  routed  to  the  sample  and  hold  circuitry  shown 
schematically  in  Rgure  A-2.  In  the  first  stage  of  this  circuit  the  sum  of  the  signals  from  the 
charge  amplifiers  is  compared  against  a  preset  voltage  threshold  level  to  determine  when  toe 
signals  are  to  be  latched.  After  toe  threshold  level  has  been  exceeded,  it  is  important  that  toe 
delay  time  of  the  triggering  signal  to  the  sample  and  hold  elements  be  small  so  that  toe 
outputs  from  toe  charge  amplifiers  will  not  saturate  before  being  latched.  Figure  9  shows  an 
oscillogram  of  toe  response  time  of  toe  nonmaskable  interrupt  (NMI)  signal.  The  lower  trace 
shows  toe  voltage  waveform  as  measured  directly  at  a  charge  collection  tab  which  denotes 
toe  time  of  charge  deposition.  The  upper  trace  shows  toe  subsequent  drop  of  toe  NMI  line. 
The  submicrosecond  response  time  of  toe  NMI  line  is  well  within  the  available  time  window  for 
sampling  the  outputs  from  toe  charge  amplifiers  as  illustrated  in  Figure  8. 

Another  important  response  time  is  toe  amount  of  time  between  toe  deposition  of  toe 
charge  and  toe  settling  of  toe  analog  division  circuitry  (Figure  A-3).  This  is,  in  fact,  a 
measurement  of  toe  total  time  required  to  make  an  analog  measurement  of  toe  charge 
deposition  position.  Figure  10  shows  toe  output  from  toe  analog  division  circuit  where  the 
voltage  waveform  from  a  collection  tab  is  again  used  to  denote  the  time  of  charge  deposition. 
The  total  time  required  from  deposition  until  the  quotient  achieves  an  asymptotic  value  is 
approximately  25  ^s.  Inasmuch  as  the  circuitry  used  here  was  constructed  from 
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generic  components  which  are  not  spedflcally  designed  for  short-cyde  times,  substantial 
reduction  in  measurement  times  should  be  possible.  The  use  of  commerdaJiy  available 
divider  modules  should  halve  the  recorded  measurement  time. 

2.3.2  Software  Evaluation.  As  illustrated  in  Figures  4  and  5.  the  uncorrected  position 
sensor  output  does  not  accurately  reproduce  the  desirable  characteristics  assodated  with  an 
idealized  position  sensor.  That  is,  the  uncorrected  output  from  the  position  sensor  is  neither  a 
highly  linear  function  of  the  actual  position  along  tfie  axis  of  determination,  nor  is  it 
independent  of  the  actual  position  along  the  orthogonal  axis.  It  is  for  these  reasons  that  the 
previously  described  correction  software  is  utilized.  In  order  to  test  the  effectiveness  of  the 
correction  software,  a  set  of  data  is  taken  for  which  20  position  measurements  are  recorded  at 
each  location  of  a  17  by  17  gnd  spanning  the  active  surface.  Figure  1 1  shows  the  mean  data 
value  at  each  position  when  the  software  correction  is  implemented.  The  correded  output 
from  the  position  sensor  shows  that  the  majority  of  the  saddle-like  distortion  associated  with 
the  uncorreded  response  of  Figure  5  has  been  eliminated.  Note  that  deviations  from  the 
idealized  inclined  plane  output  charaderistics  are  predominantly  localized  to  the  perimeter  of 
the  active  surface. 

Figure  12  simultaneously  illustrates  the  correded  response  from  both  axes  of  the  position 
sensor.  Comparison  with  the  uncorrected  response,  of  Figure  7  shows  that  the  data 
correction  software  has  significantly  improved  the  correlation  between  the  measured  and  ideal 
response  as  illustrated  in  Figure  6. 

In  order  to  quantitatively  discuss  the  amount  of  distortion  associated  with  the  uncorreded 
and  correded  position  sensor  response,  it  is  convenient  to  consider  the  residual  values 
associated  with  each  case.  In  this  context,  the  residual  value  is  defined  to  be  the  difference 
between  the  actual  location  of  charge  deposition  and  the  mean  location  as  determined  by  the 
position  sensor.  A  residual  value  of  zero  therefore  indicates  an  exad  determination  of 
position  with  the  magnitude  of  nonzero  residuals  reflecting  the  amount  of  error,  or  lack  of 
accuracy,  associated  with  the  measurement. 

Figure  13  illustrates  the  residuals  associated  with  the  uncorreded  position  sensor  data  in 
two  different  formats.  On  the  top,  a  three-dimensional  plot  of  residuals  graphically  displays 
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Figure  13.  Residual  Values  for  Uncorrected  X  Axis  Position  Measurements  With 
Three-Dimensional  Representation  at  the  Too  and  Two-Dimensional 

Representation  at  the  Bottom. 
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Figure  14.  Residual  Values  For  X  Axis  Software-Corrected  Position  Measurements  With 
Three-Dimensional  Representation  at  the  Top  and  Two-Dimensional 
Representation  at  the  Bottom. 
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the  relative  magnitude  and  sense  of  the  errors  associated  with  each  measurement  location. 
This  type  of  plot  is  useful  for  conveying  an  overall  image  of  where  the  accuracy  of  the  system 
is  limited.  In  this  type  of  representation,  an  idealized  position  sensor  has  a  three-dimensional 
residual  plot  which  is  a  plane  with  an  X  residual  value  of  zero.  For  a  more  quantitative 
assessment,  a  two-dimensional  plot  of  the  same  data  is  more  useful,  as  shown  on  the  lower 
figure.  In  this  illustration,  the  actual  X  position  axis  has  been  compressed  by  plotting  the  X 
residual  value  vs.  the  actual  Y  position  for  the  family  of  actual  X  positions.  The  figure  on  the 
bottom  can  also  be  thought  of  as  the  figure  on  the  top  when  viewed  from  a  vantage  point 
located  a  iarge  distance  down  the  actual  X  position  axis.  The  envelope  defined  by  this  family 
of  c  rves  represents  the  "worst-case"  errors  associated  with  the  uncorrected  X  position 
measurements.  X,.  at  each  actual  Y  position.  This  two-dimensional  plot  clearly  shows  that  the 
greatest  errors  in  the  uncorrected  X  position  measurement  occur  for  minimum  and  maximum 
values  of  the  actual  Y  position.  As  previously  discussed,  the  enhanced  error  along  the  edges 
of  the  active  surface  is  due  to  the  presence  of  the  charge  collection  tabs  used  for 
simultaneous  measurements  of  the  position  in  the  orthogonal  direction.  The  worst-case  error 
for  the  uncorrected  X  axis  data  is  located  in  this  region  with  a  magnitude  of  9.2  parts  out  of  40 
or  23%  of  full  scale.  Combining  the  X  and  Y  axis  data,  the  worst-case  error  for  the 
uncorrected  X,,Y,  data  is  12.0  parts  out  of  40  or  30%  of  full  scale. 

The  effective  accuracy  of  the  position  sensor  can  be  improved  by  limiting  the  functional 
region  of  the  active  surface  to  an  interior  area.  This  accuracy  increase  is  obviously 
associated  with  a  decrease  in  position  sensing  coverage.  If  the  working  active  surface  is 
taken  to  be  the  area  defined  by  the  inner  15  by  15  grid  of  data  points,  instead  of  the  entire  17 
by  17  grid  of  the  active  surface,  the  usable  sensing  area  is  reduced  to  76%.  In  this  interior 
region,  the  worst-case  error  for  the  uncorrected  X  axis  data  drops  to  6.7  parts  out  of  40  or 
17%  of  full  scale.  The  worst-case  error  for  the  combined  X,,Y,  data  is  reduced  to  7.9  parts  out 
of  40  or  20%  of  full  scale. 

The  two-dimensional  plot  of  Figure  13  highlights  the  worst-case  data  points  but  offers  little 
quantitative  information  about  the  average  error  associated  with  the  position  measurements. 
Numerical  analysis  of  the  full  area  sensor  data  yields  an  average  residual  magnitude  for  the 
uncorrected  X  axis  data  of  2.5  parts  out  of  40  or  6.3%  of  full  scale.  For  the  combined  X„Y, 
data,  the  average  residual  magnitude  is  4.0  parts  out  of  40  or  10%  of  full  scale.  The  average 
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residual  magnitude  for  the  X  axis  data  in  the  previously  defined  interior  region  drops  to  2.0 
parts  out  of  40  or  5%  of  full  scale.  In  this  same  region,  the  average  residual  magnitude  of  the 
combined  X^Y,  data  is  3.2  parts  out  of  40  or  8%  of  full  scale. 

Similar  plots  of  the  residual  values  for  software-corrected  measurements  are  shown  in 
Figure  14.  The  three-dimensional  plot  on  the  top  can  be  seen  to  more  closely  approximate 
the  planar  residual  plot  of  an  idealized  position  sensor.  On  the  bottom,  the  reduced  envelope 
of  the  two-dimensional  plot  shows  that  the  worst-case  accuracy  errors  have  been  significantly 
reduced.  For  the  full  area  sensor,  the  worst-case  error  for  the  comected  X  axis  data  is  5.1 
parts  out  of  40  or  13%  of  full  scale.  The  average  residual  magnitude  for  the  X  axis  data  over 
this  full  area  is  0.55  parts  out  of  40  or  1 .4%  of  full  scale.  The  combined  X^Y^  data  shows  a 
worst-case  error  of  7.1  parts  out  of  40  or  18%  of  full  scale  and  an  average  residual  magnitude 
of  0.9  parts  out  of  40  or  2.3%  of  full  scale.  Once  again,  the  most  significant  accuracy  errors 
are  localized  along  the  perimeter  of  the  active  surface.  Therefore,  the  effective  accuracy  of 
the  position  sensor  can  again  be  increased  at  the  expense  of  usable  sensing  area.  For  the 
previously  defined  interior  area,  the  worst-case  error  for  the  corrected  X  axis  data  drops  to  2.0 
parts  out  of  40  or  5%  of  full  scale  while  the  average  residual  magnitude  is  reduced  to  0.36 
parts  out  of  40  or  0.9%  of  full  scale.  In  this  same  region,  the  combined  data  worst-case 
error  drops  to  3.6  parts  out  of  40  or  9%  of  full  scale  while  the  average  residual  magnitude  is 
reduced  to  0.63  parts  out  of  40  or  1 .6%. 

The  utility  of  ttie  software  correction  process  in  reducing  accuracy  errors  is  revealed  by 
comparing  the  errors  for  the  uncorrected  and  corrected  cases.  For  the  full  area  sensor,  the 
correction  routines  reduce  the  absolute  worst-case  error  for  the  X  axis  data  from  23%  to  13% 
and  the  average  residual  magnitude  from  6.3%  to  1.4%.  Over  the  full  area  of  the  sensor,  the 
correction  routines  reduce  the  combined  X  and  Y  data  worst-case  error  from  30%  to  1 8%  and 
the  average  residual  magnitude  from  10%  to  2.3%.  If  only  the  interior  76%  of  the  active 
surface  area  is  utilized,  the  software  correction  drops  the  X  axis  data  worst-case  error  from 
1 7%  to  5%  while  the  average  residual  magnitude  decreases  from  5%  to  0.9%.  The  effect  of 
the  software  correction  on  the  combined  X  and  Y  axis  data  in  the  interior  region  is  to  reduce 
the  worst-case  error  from  20%  to  9%  and  the  average  residual  magnitude  from  8%  to  1 .6%. 
Thus,  the  correction  routines  decrease  the  worst-case  error  by  approximately  a  factor  of  two 
and  the  average  residual  magnitudes  by  a  factor  of  four  or  five. 
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The  ptBvious  discussion  and  quantitative  anaiysis  of  the  calibration  of  the  position  sensor 
utilizes  calibration  and  test  data  comprised  of  the  mean  values  of  an  ensemble  of 
measurements  at  each  tested  location.  Therefore,  these  determined  accuracies  are  an 
indication  of  the  average  error  between  the  actual  and  measured  positions.  In  addition  to  the 
accuracy  associated  with  position  measurements,  it  is  also  important  to  consider  the  precision, 
or  reproducibility,  of  these  measurements.  This  precision  measurement  describes  the  range 
of  values  over  which  the  ensemble  of  measurements  is  spread.  Rgure  15  qualitatively 
displays  the  reproducibility  of  five  sets  of  X  value  data.  These  data  were  collected  by 
recording  500  sequential  raw  X  values  at  5  different  locations  on  the  position  sensor's  active 
surface.  The  data  in  sets  A  B.  C,  D.  and  E  correspond  to  actual  X  positions  of  7. 14,  21 , 28. 
and  34.  A  centrally  located  actual  Y  position  was  used  for  all  of  these  tests.  In  general,  the 
distributions  for  these  measurements  are  not  symmetrical.  This  is  particularly  true  for  mea¬ 
surements  taken  near  an  edge  of  the  active  surface.  Therefore,  these  data  are  not  nomnaily 
distributed,  and  a  simple  standard  deviation  should  not  be  used  to  represent  the  range  of  data 
values.  Instead,  the  precision  is  quantified  by  ranking  the  individual  measurements  and 
specifying  the  range  of  values  over  which  a  given  percentage  of  the  data  is  contained.  For 
these  test  data,  the  range  of  values  which  contains  90%  of  the  individual  measurements  will 
be  considered.  Since  20  data  values  were  recorded  at  each  location,  this  is  equivalent  to 
removing  the  lowest  and  highest  values  and  determining  the  span  in  position  of  the  18 
centrally  ranked  measurements. 

The  magnitudes  of  the  spans  of  X  axis  values  which  contain  90%  of  the  software- 
corrected  measurements  are  illustrated  graphically  in  Rgure  1 6.  The  formats  of  the  plots  and 
the  magnitudes  of  the  axes  scales  are  similar  to  previous  plots  to  facilitate  direct  comparison. 
Once  again,  the  largest  errors  are  localized  along  the  outer  edges  of  the  active  surface.  The 
largest  recorded  span  in  X  values  which  contained  90%  of  the  individual  measurements  over 
the  full  area  sensor  was  12  parts  out  of  40  or  30%  of  full  scale.  However,  the  average  X 
value  span  for  the  full  area  sensor  was  much  lower  at  1 .74  parts  out  of  40  or  4.4%  of  full 
scale.  Over  the  previously  defined  interior  area,  the  largest  X  value  span  was  8  parts  out  of 
40  or  20%  of  full  scale  with  an  average  value  of  1 .44  parts  out  of  40  or  3.6%  of  full  scale. 

The  combined  deviations  of  the  individual  position  measurements  along  the  X  and  Y  axes 
are  illustrated  in  Figure  17.  For  each  position  of  the  17  by  17  grid  of  test  locations,  the  area 
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Rgure  15.  Histoqran  of  Raw  X  Values. 

of  the  position  sensor  which  contained  90%  of  the  individual  measurements  is  outlined.  Small 
rectangles  indicate  a  close  grouping  of  individual  data  points  while  larger  rectangles  mark 
areas  of  greater  deviation.  Once  again,  this  illush’ation  highlights  the  fact  that  the  largest 
deviations  are  located  along  the  edges  of  the  transducer  surface. 

Depending  on  the  sense  of  a  deviation  with  respect  to  the  sense  of  the  associated 
residual,  a  deviation  from  the  mean  position  can  either  move  a  recorded  position  closer  to,  or 
farther  away  from,  the  actual  location.  To  consider  the  worst-case  situation,  it  is  assumed  that 
the  sense  of  the  deviation  is  the  same  as  the  sense  of  the  residual  so  that  these  two  sources 
of  error  will  add.  Thus,  a  conservative  measure  of  the  overall  performance  that  can  be 
expected  from  the  position  sensor  is  to  consider  the  distance  between  the  actual  ^arge 
deposition  location  and  the  most  distant  comer  of  the  associated  deviation  rectangle,  as 
displayed  in  Figure  17.  This  distance  is  termed  the  "position  error."  Using  this  distance  to 
define  the  radius  of  a  drcle,  at  least  90%  of  all  individual  position  measurements  can  be 
expected  to  fall  within  this  circle  when  it  is  centered  at  the  actual  charge  deposition  location. 
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Positions  With  Three-Dimensional  Reoresentatation  at  the  Tod  and  Two- 
Dimensional  Representation  at  the  Bottom. 
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Figure  18  plots  the  confidence  level  that  90%  of  all  randomly  located  individual  position 
measurements  will  be  within  a  given  position  error  of  the  actual  charge  deposition  location. 
Data  are  presented  for  both  the  full  and  reduced  sensing  areas.  This  plot  was  constructed 
from  the  ranked  position  errors  of  the  289  positions  of  the  17  by  17  grid  of  test  locations.  The 
ranks  were  normalized  from  0%  to  100%  to  form  the  abscissa  while  the  ranked  position  errors 
formed  the  ordinate.  For  example,  the  largest  observed  position  error  for  the  reduced  area 
sensor  was  14.5%  of  full  scale.  Therefore,  the  confidence  level  is  100%  that  randomly 
located  position  measurements  will  yield  individual  position  measurements  such  that  90%  will 
be  wittiin  14.5%  of  full  scale  of  the  actual  position.  Similarly,  the  confidence  level  is  80%  that 
randomly  located  position  measurements  in  the  reduced  area  will  have  90%  of  the  individual 
measurements  within  6.3%  of  fuil  scale  of  tire  actual  position.  Selected  confidence  levels  and 
the  associated  position  errors  are  tabulated  in  Table  1. 

2.3.3  Simulated  Reid  Testing.  For  laboratory  calibration  and  perfonnance  evaluation, 
a  metallic  cover  plate  was  used  which  had  an  array  of  small  holes  spanning  the  active  surface 
of  the  position  sensor.  Charge  could  be  repeatedly  deposited  at  well-defined  locations  by 
passing  a  fine  metallic  probe  through  these  holes  and  shorting  the  cover  plate  to  the  active 
surface.  For  field  applications,  the  stored  charge  would  be  deposited  in  a  somewhat  more 
violent  and  less  controlled  manner.  The  impacting  object  would  deforni  the  cover  plate  at  the 
point  of  impact  where  the  stored  charge  would  be  transferred  to  the  active  surface.  These  two 
charge  deposition  methods  are  expected  to  be  equivalent  because,  as  previously  stated,  it  is 
the  ratio  of  charges  collected  at  the  various  tabs  which  is  critical  and  not  the  magnitude  of  the 
deposited  charge.  To  test  the  sensitivity  of  position  determination  to  the  actual  amount  of 
deposited  charge,  a  series  of  measurements  was  taken  using  charge  storage  capacitors  of 
different  magnitudes.  Position  measurements  were  recorded  at  five  different  locations  while 
the  charge  storage  capacitance  was  varied  over  six  orders  of  magnitude.  Twenty  individual 
position  measurements  were  logged  at  each  position  and  capacitance  value.  The  X  and  Y 
mean  positions,  5%  rank  positions,  and  95%  rank  positions  are  listed  in  Table  2  for  the 
various  capacitance  values.  As  expected,  neither  the  accuracy  nor  the  reproducibility  of  the 
position  measurements  appear  to  be  affected  by  the  variation  in  the  magnitude  of  the 
deposited  charge. 
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Rgure  18.  Position  Error  vs.  Confidence  Level. 


Table  1 .  Dependence  of  Confidence  Level  on  Position  Error 


Confidence  Level 
(%) 

Full  Error  Position  Level 
(%  of  FS) 

Reduced  Area  Position 
Error 
(%  of  FS) 

10 

1.9 

1.8 

20 

2.5 

2.3 

30 

3.1 

2.8 

40 

3.9 

3.4 

50 

4.5 

4.0 

60 

5.3 

4.5 

70 

6.3 

5.1 

80 

7.4 

6.3 

90 

11.2 

7.3 

100 

23.5 

14.5 

Of  course  the  best  way  to  evaluate  the  field  perfonnance  of  the  position  detector  is  to  fire 
projectiles  at  it.  This  allows  the  system  to  be  tested  when  charge  is  transferred  to  the  active 
surface  following  the  deformation  of  the  cover  plate.  Projectiles  of  4.5*mm  diameter  were 
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impacted  on  the  position  sensor  with  a  nomined  velocity  of  100  m/s.  The  force  of  these 
impacts  was  sufficient  to  deform  and  pierce  the  cover  plate  with  the  associated  deposition  of 
stored  charge.  Table  3  lists  the  recorded  analog  signal  measurement  times;  raw.  corrected, 
and  actual  X  and  Y  positions,  and  position  errors.  The  impact  sites  were  spread  over  the 
position  sensor  active  surface  in  a  roughly  homogenous  pattern.  Shots  with  impact  locations 
outside  of  the  previously  defined  interior  region  are  denoted  by  an  asterisk  next  to  the  test 
number.  Ten  of  the  37  test  shots,  or  27%,  fell  outside  this  interior  region.  Since  the  external 
region  comprises  24%  of  the  active  surface,  this  simulated  field  test  represents  a  survey  of  the 
entire  active  surfoce.  Figure  19  compares  the  full  area  confidence  level  statistics  as 
determined  in  the  laboratory  against  the  field  test  results.  The  laboratory  full  area  position 
errors  include  only  90%  of  the  individual  position  measurements  whereas  the  single-shot 
statistics  of  the  field  tests  effectively  represent  100%  of  each  shot.  In  spite  of  the  more 
stringent  definition  for  the  field  test  position  error,  the  confidence  level  for  a  given  position 
error  is  higher  for  the  field  test  data  than  the  laboratory  simulation.  Also,  the  overall  shapes  of 
these  two  curves  are  quite  similar.  This  confirms  the  previous  assertion  that  the  performance 
statistics  developed  in  the  laboratory  are  a  conservative  estimate  of  the  position  sensors  cap¬ 
abilities.  Thus,  the  performance  of  the  position  sensor  as  measured  in  the  laboratory  is  a 
valid  estimate  of  the  field  performance  of  this  device. 

Another  important  perfonnance  parameter  for  an  impact  position  sensor  is  the  ’multiple-hit” 
capability.  That  is,  how  is  the  performance  of  the  sensor  degraded  by  the  damage  sustained 
during  previous  impacts?  For  the  previously  described  field  testing,  a  thick  hard  back  plate 
was  located  directly  behind  the  electrically  resistive  surface.  Thus,  although  the  resistive 
surface  was  damaged  by  each  test,  the  active  surface  was  not  actually  pierced.  To  study  the 
effect  of  active  surface  deterioration  on  position  sensing,  a  hole  was  cut  into  the  active 
surface,  and  the  laboratory  performance  evaluation  was  repeated.  A  12-mm-diameter  hole 
was  cut  into  the  resistive  surface  which  yields  a  hole  size  that  is  6%  of  the  active  surface 
edge  length.  This  hole  was  positioned  so  that  it  did  not  coincide  with  any  of  the  positions  of 
the  17  by  17  grid  of  test  sites.  Figure  20  illustrates  the  relative  positions  of  the  test  sites  and 
the  1 2-mm  hole.  Superimposed  on  this  figure  are  the  corresponding  corrected  positions. 
Comparison  of  Figure  20  with  the  corrected  positions  for  the  undamaged  sensor,  shown  in 
Figure  12,  reveals  that  the  position  sensor  output  characteristics  have  changed  very  little. 
Quantitative  analysis  confirms  this  observation.  The  introduction  of  the  hole  in  the  active 
surface  increases  the  full  area  average  residual  magnitude  from  2.3%  to  2.6%  while  the 
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Rgure  19.  Position  Error  vs.  Confidence  Levet  for  Simulated  Reid  Test  Data. 


Figure  20.  Corrected  Position  Sensor  X  and  Y  Output  Characteristics  for  a  Damaged 
Active  Surface. 
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reduced  area  average  residual  magnitude  changes  from  1.6%  to  2.3%.  These  small 
increases  in  average  residual  magnitude  indicate  that  this  impact  position  sensor  exhibits 
potential  for  substantial  muitihit  capabilities. 

2.4  One-Dimensional  Measurements. 

2.4.1  Position  Sensing  Ti-ansducer.  One  source  of  the  nonlinearity  in  the  raw  response  of 
this  impact  sensor  transducer,  and  the  subsequent  need  for  data  correction,  is  the  presence  of 
the  charge  collection  tabs  which  are  simultaneously  measuring  the  position  in  the  orthogonal 
direction.  A  logical  extension  of  the  results  already  presented  is  to  consider  the  performance 
of  the  impact  sensor  with  only  a  single  pair  of  charge  collection  tabs,  as  illustrated  in 

Figure  21.  For  this  configuration,  the  remaining  two  charge  collection  tabs  are  extended  to 
span  the  entire  edge  of  the  active  surface  to  reduce  edge  effects  and  more  closely 
approximate  an  infinitely  long  one-dimensional  system.  Of  course,  this  type  of  arrangement 
can  only  provide  position  information  along  tiie  axis  which  is  normal  to  the  collection  tabs. 
However,  a  second  transducer  surface  with  the  collection  tabs  oriented  in  the  orthogonal 
direction  could  be  used  to  provide  impact  information  along  the  second  dimension  for  a 
complete  two-dimensional  position  determination.  The  added  complication  of  the  second 
transducer  surface  may  be  more  than  offset  by  the  potential  for  projectile  velocity 
measurements.  By  measuring  the  time  interval  between  the  impacts  of  the  projectile  on  the 
two  transducers,  and  knowing  the  separation  distance  between  the  transducers,  the  projectile 
velocity  can  be  determined. 

2.4.2  One-Dimensional  Preliminary  Results.  Figure  22  illustrates  the  raw  output 
characteristics  for  the  one-dimensional  position  sensor.  In  contrast  to  the  raw  output 
characteristics  of  the  transducer  surface,  which  simultaneously  measured  positions  in  two 
dimensions  (shown  in  Figure  5),  the  one-dimensionai  transducer  response  is  quite  linear.  The 
raw  output  from  the  one-dimensional  transducer  does  not  completely  span  the  range  of 
determined  values.  One  method  for  correcting  this  deficiency  is  to  employ  the  previously 
described  software  correction  techniques.  The  one-dimensional  corrected  output 
characteristics  are  illustrated  in  Figure  23.  However,  the  true  worth  of  the  one-dimensional 
measurement  scheme  may  lie  in  the  potential  for  completely  eliminating  the  need  for 
computer-controlled  data  correction.  By  modifying  the  output  stage  electronics  to  expand  the 
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Figure  22.  Uncorrected  One-Dimensional  Position  Sensor  Output  Characteristics. 


Rgure  23.  Corrected  One-Dimensional  Position  Sensor  Output  Characteristics. 

available  zero-level  and  amplification  control,  the  linear  raw  output  from  the  one-dimensionai 
transducer  can  be  spread  out  over  the  available  range  of  determined  values.  Elimination  of 
the  need  for  software  data  correction  will  increase  the  speed  of  the  position  measurement 
while  decreasing  the  cost  and  complexity  of  the  entire  system. 

3.  SUMMARY 

A  proof-of-principle  integrated  hardware-software  system  has  been  developed  which 
rapidly  measures  impact  locations  on  a  planar  surface.  The  transducer  for  this  position 
sensor  consists  of  an  electrically  resistive  square  sheet  with  charge  collection  tabs  along  the 
outer  edges.  When  electrical  charge  is  deposited  at  the  point  of  impact  the  rate  of  charge 
transfer  to  each  collection  tab  depends  on  the  resistance  and,  therefore,  the  physical  distance 
between  the  point  of  charge  deposition  and  the  tab.  The  impact  position  is  deduced  from  the 
rate  at  which  charge  is  transferred  to  the  collection  tabs.  Two  transducer  geometries  have 
been  investigated.  Predominantly,  a  four-tab  arrangement  has  been  tested  which  allows  a 
single  transducer  to  simultaneously  determine  position  along  two  axes.  In  addition,  a  cursory 
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investigation  has  considered  a  transducer  geometry  with  two  charge  coilection  tabs,  which 
allows  position  to  be  measured  along  a  single  axis. 

The  electronic  hardware  which  is  used  to  process  the  collected  charge  has  been  tested 
and  is  documented  in  Appendix  A.  To  increase  the  accuracy  of  position  measurements,  a 
collection  of  software  routines  which  are  used  to  digitize  the  output  signal,  precalibrate  the 
system,  and  correct  the  raw  position  measurements  has  been  developed.  These  software 
routines  are  described  and  listed  in  Appendix  B.  Application  of  the  software  correction  utilities 
is  shown  to  reduce  accuracy  errors  by  a  factor  of  four.  Relevant  accuracy  and  precision 
parameters  are  tabulated  below  for  errors  along  a  single  axis.  Accuracy  errors  are  also 
tabulated  for  combined  two  axis  measurements.  This  data  was  obtained  using  the  four-tab 
transducer  geometry  in  a  laboratory  environment.  In  Table  4.  *full  area”  indicates  that  the 
entire  active  surface  of  the  position  sensor  was  utilized  while  ‘reduced  area*  refers  to  tests  in 
which  only  the  central  76%  of  the  sensor’s  active  area  was  tested.  Accuracy  error  refers  to 
the  difference  between  the  actual  position  and  the  mean  position  of  a  large  number  of 
individual  measurements.  Precision  error  is  a  measure  of  the  range  of  values  over  which  the 
ensemble  of  individual  measurements  is  spresuf.  For  these  tests,  the  precision  error  was 
quantified  by  considering  the  range  of  values  which  contained  90%  of  the  individual 
measurements.  All  error  values  are  given  as  a  percentage  of  the  active  surface  edge  length. 

For  the  combined  two-dimensional  position  measurements,  the  notions  of  accuracy  and 
precision  error  were  combined  into  a  single  gauge  of  total  error.  Empirical  laboratory  tests 
determined  the  confidence  levels  that  90%  of  all  individual  position  measurements  would  fail 
within  a  circle  of  a  given  radius  or  position  error.  Simulated  field  testing  of  the  position  sensor 
yielded  total  errors  which  were  equal  to  or  smaller  than  those  predicted  by  laboratory 
evaluation.  The  apparent  equivalence  between  the  performance  of  the  position  sensor  in 
laboratory  and  field  environments  is  due  in  part  to  the  insensitivity  of  the  system  to  the  actual 
amount  of  deposited  charge.  This  was  demonstrated  by  varying  the  capacity  of  the  charge 
storage  capacitor  over  six  orders  of  magnitude  with  no  variation  in  system  performance.  The 
ability  to  predict  field  performance  through  laboratory  testing  is  a  particularly  convenient 
capability. 
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Table  4.  Position  Sensor  Error  Parameters 


1-Axis 

Maximum 

Accuracy 

Error 

(%) 

1-Axis 

Average 

Accuracy 

Error 

(%) 

1-Axis 

Maximum 

Precision 

Error 

(%) 

I-Axis 

Average 

Precision 

Error 

(%) 

2-AxiS 

Maximum 

Accuracy 

Error 

(%) 

2-Axis 

Average 

Accuracy 

Error 

(%) 

Uncorrect 
Full  Area 

23 

6.3 

— 

— 

30 

10  1 

Uncorrect 

Reduced 

Area 

17 

5.0 

— 

20 

8.0  1 

Correct 

Full  Area 

13 

1.4 

30 

4.4 

18 

2.3 

Correct 

Reduced 

Area 

m 

0.9 

20 

3.6 

9.0 

1.6 

Laboratory  evaluation  was  also  utilized  in  a  study  of  the  system’s  ability  to  perform  after 
sustaining  transducer  damage.  A  hole  was  cut  into  the  active  surface,  and  the  error 
parameters  were  remeasured.  The  system  displayed  a  minimal  reduction  in  perfomiance 
which  indicates  that  for  moderate  damage,  the  impact  position  sensor  will  display  considerable 
multihit  capabilities. 

The  preliminary  investigation  of  a  one-dimensional,  two-charge  collection  tab  geometry 
revealed  a  conveniently  linear  transducer  response.  With  improved  signal  processing,  this 
system  may  obviate  the  need  for  complicated  data  correction.  In  order  to  obtain  two- 
dimensional  position  information,  two  of  these  one-dimensional  transducers  would  be  required. 
The  addition  of  a  second  transducer  surface  offers  the  potential  for  additional  projectile 
velocity  measurements. 

It  is  stressed  that  the  impact  position  sensor  described  in  this  report  is  a  proof-of-principle 
device  and,  therefore,  does  not  exemplify  the  ultimate  performance  which  might  be  achieved 
by  a  refined  device.  In  fact,  it  is  clear  that  the  performance  of  this  position  detector  can  be 
improved  in  a  number  of  ways.  Modeling  of  the  charge  transfer  process  on  the  active  surface 
may  suggest  optimizing  configurations  for  the  resistive  sheet  and  collection  tabs  to  improve 
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the  systems  accuracy.  Refinements  in  the  signal  processing  electronics  may  reduce  the  error 
associated  with  measurement  reproducibility.  The  time  required  to  make  the  position 
measurement  can  be  reduced  in  a  number  of  ways.  Analog  and  digital  signal  response  times 
can  both  be  halved  by  utilizing  high-speed  electronic  components.  And.  finally,  the  time 
required  for  the  data  correction  process  can  be  reduced  by  employing  correction  look-up 
tables  which  are  stored  in  hardwired  programmable  read  only  memory  (PROM)  chips  rather 
than  software-controlled  computer  iook-up  tables. 
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APPENDIX  A: 

SIGNAL  PROCESSING  CIRCUITRY 
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A.1  Charofl  Amolifief.  A  schematic  for  the  dual-channel  charge  amplifiers  is  shown  in 
Figure  A-1.  Charge  from  opposing  tabs  on  the  position-sensing  surtece  is  feed  into  IN,  and 
IN,  to  the  inverting  inputs  of  the  operational  amplifier.  Since  the  noninverting  inputs  are  tied  to 
ground,  the  op  amp  outputs  will  try  to  maintain  the  inverting  inputs  at  ground  potential  by 
feeding  charge  of  the  opposite  sense  through  the  selected  feedback  capacitor  to  the  inverting 
input  As  charge  passes  through  the  feedback  capacitor,  a  potential  develops  across  the 
capacitor  of  magnitude. 


AV^»AQ/C.  (A-1) 

where  V  is  the  potential  developed  across  capacitance  C  for  an  amount  of  deposited  charge 
Q.  From  this  equation,  it  can  be  seen  that  the  magnitude  of  the  output  voltage  can  be  scaled 
by  the  magnitude  of  the  feedback  capacitor.  Thus,  the  amplification  of  these  charge  amplifiers 
is  controlled  by  selection  of  a  suitable  feedback  capacitor. 

A  large  value  resistor  is  included  in  parallel  across  the  feedback  capacitors  to  prevent  the 
effective  integration  of  the  op  amp's  input  bias  current  which,  although  small,  would  eventually 
saturate  the  amplifier.  The  use  of  a  feedback  resistor  causes  the  output  of  the  amplifier  to 
"droop"  with  time.  However,  for  the  values  of  capacitors  and  resistors  shown,  the  hold-time  of 
the  charge  amplifier  is  sufficiently  long  to  allow  the  output  to  be  captured  by  the  subsequent 
sample  and  hold  circuitry.  In  addition,  the  relative  zero  levels  of  the  two  channels  can  be 
adjusted  by  the  potentiometer  connected  to  the  op  amp  balance  pins. 

A.2  Sample  and  Hoid  Circuitry.  A  schematic  of  the  sample  and  hold  is  shown  in 
Figure  A-2.  The  output  signals  from  the  charge  amplifiers  are  fed  into  this  circuit  as  SIGNAL 
1  and  SIGNAL  2.  These  input  signals  are  processed  by  the  circuitry  of  IC3  to  yield  a  voltage 
which  is  the  sum  of  the  two  input  signals.  This  summation  voltage  is  compared  against  a 
preset  threshold  voltage  at  IC4.  When  the  summation  voltage  exceeds  the  threshold  voltage, 
the  one-shot  of  ICS  is  triggered.  This  trigger  pulse  serves  a  number  of  purposes  and  is 
therefore  buffered  using  transistor  T1.  One  task  of  the  trigger  signal  is  to  illuminate  a  light 
emitting  diode  which  notifies  the  operator  that  an  event  has  occurred.  The  trigger  signal  also 
pulls  down  an  active-low  nonmaskable  interrupt  (NMI)  line  which  can  be  used  to  synchronize 
associated  computer-controlled  data  acquisition  equipment.  The  computer  used  for  this  task 
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operates  with  0-V  to  5-V  TTL  signal  levels.  However,  this  sample  and  hold  circuitry  operates 
with  +V  and  -V  voltages  of  +12  and  -12  V.  Therefore,  the  trigger  pulse  is  passed  through  a 
voltage  divider  for  use  as  an  NMI  signal.  Rnaiiy,  the  trigger  pulse  is  also  fed  to  IC6  and  IC7 
which  sample  and  hold  the  signals  from  the  charge  amplifiers.  The  captured  signal  from 
SIGNAL  2  is  output  as  SIGNAL  2H.  In  addition,  SIGNAL  2H  is  added  to  the  captured  signal 
SIGNAL  1H  by  ICS  and  output  as  SIGNAL  2H+SIGNAL  1H. 

A.3  Analog  Division  Circuit  The  division  of  SIGNAL  2H  by  SIGNAL  2H+SIGNAL  1H  is 
performed  by  taking  the  logarithms  of  these  voltages,  subtracting  the  logarithms,  and  then 
determining  the  antilogarithm  of  the  subtraction.  Such  an  approach  is  made  possible  by  the 
convenient  fact  that  for  a  silicon  transistor,  the  magnitude  of  the  base-emitter  voltage,  is  a 
logarithmic  function  of  the  collector  current,  l^,  ttirough  the  transistor.  The  circuitry  which 
performs  these  tasks  is  shown  schematically  in  Figure  A-3  (Jung  1980). 

This  circuit  operates  as  follows.  Since  the  noninverting  inputs  of  op  amps  IC9  and  IC10 
are  tied  to  ground,  the  associated  inverting  inputs  are  maintained  at  an  effective  ground 
potential.  Therefore,  SIGNAL  2H  sets  up  a  current  flow  through  its  input  resistor  which  is 
equal  to  (SIGNAL  2H)/R.  This  input  current  is  equal  to  the  collector  current  through  transistor 
Tj.  By  the  logarithmic  relationship  between  I,,  and  it  follows  that  the  voltage  drop  between 
the  base  and  emitter  of  Tj  must  have  a  magnitude  which  is  proportional  to  ln[(SIGNAL  2H)/R]. 
Because  Tj  is  an  npn  transistor  and  the  base  is  tied  to  ground,  it  follows  that  the  sense  of  the 
voltage  at  the  emitter  of  this  transistor  will  be  negative  with  respect  to  ground.  By  similar 
reasoning,  the  current  through  transistor  T3  is  (SIGNAL  2H+SIGNAL  which  generates  a 
voltage  between  the  base  and  emitter  of  ln[(SIGNAL  2H+SIGNAL  1H)/RJ.  The  emitter  of  T3  is 
tied  to  the  emitter  of  Tj  which  implies  that  the  magnitude  of  the  potential  at  the  emitter  of  T3 
must  be  ln(SIGNAL  2H/R).  Since  T3  is  an  npn  transistor,  the  conditions  that  the  emitter  is  at  a 
potential  of  magnitude  ln(SIGNAL  2H/R),  and  the  voltage  drop  between  the  base  and  emitter 
has  a  magnitude  of  ln[(SIGNAL  2H+SIGNAL  1H)/R]  requires  that  the  base  be  at  a  potential  of 
ln(SIGNAL  2H/R)  -  ln[(SIGNAL  2H+SIGNAL  1H)/RJ.  This  base  potential  is  provided  by  output 
pin  7  of  IC9  and  represents  the  difference  between  the  logarithms  of  the  two  input  signals.  If 
the  antilogarithm  of  this  signal  were  taken,  the  required  division  would  be  completed. 

However,  for  convenience  this  signal  is  combined  with  another  adjustable  signal,  referred  to 
as  SCALE,  which  is  on  the  order  of  10  V  and  can  be  used  to  vary  the  scale  of  the  output 
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Rgure  A-3.  Schematic  of  Analog  Division  Circuit 


signal.  Transistor  is  configured  much  like  T,.  Since  the  collector  current  through  is 
SCALE/R,  and  the  base  potential  is  set  at  ln(SIGNAL  2H/R)  -  ln[(SlGNAL  2H-^SIGNAL 
it  follows  that  the  emitter  must  be  at  a  potential  of  magnitude  ln(SCALE/R)  +  ln(SIGNAL  2H/R) 
-  ln[(SIGNAL  2H4SIGNAL  1H)/R].  Using  the  laws  of  logarithms  and  the  fact  that  all  the  fls  are 
equal,  this  can  be  rewritten  as: 


SIGNAL  2H 

SIGNAL  2H-*-SIGNAL  ' 

Thus,  all  that  is  left  to  be  done  is  to  take  the  antilogarithm  of  the  above  expression,  and  the 
necessary  division  will  have  been  performed,  complete  with  a  variable  scaling  factor.  This 
task  is  accomplished  by  transistor  T,.  Because  the  base  of  this  transistor  is  grounded,  the 
base  to  emitter  voltage  drop  will  have  a  magnitude  equal  to  the  above  expression.  A  of 
this  magnitude  allows  a  collector  current  of  magnitude  equal  to  the  antilogahthm  of  which 
is  SCALE  •  SIGNAL  2H/(SIGNAL  2H+SIGNAL  1H).  This  current  passes  through  the  feedback 
resistor  of  the  associated  op  amp  which  has  a  value  which  is  only  one-tenth  that  of  the  input 
resistors.  Therefore,  the  output  signal  is  reduced  by  a  factor  of  ten,  which  offsets  the 
approximate  factor  of  ten  introduced  by  the  SCALE  signal,  and  the  analog  division  is 
complete. 


In  \  SCALE  * 
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The  various  feedback  capacitors  are  included  for  stability  against  osditation  while  the 
damp  diodes  protect  the  transistors  against  reverse  biasing.  For  maximum  temperature 
stability,  the  transistor  pairs  should  be  either  monolithic  transistor  pairs,  such  as  LM394’s  or 
well-matched  discretes. 

A.4  Voltaae  Amplifier.  A  differential  input  voltage  amplifier  is  induded  as  the  iast  step  of 
analog  signal  processing.  The  abiiity  to  vary  both  the  DC  offset  and  the  signal  amplitude 
allows  the  final  output  signal  to  be  tailored  for  further  digital  signal  processing  or  dispiay. 
Through  the  use  of  differential  amplifiers,  common-mode  noises,  such  as  a  60-Hz  "hum"  and 
capadtiveiy  coupled  microprocessor  docking  signals,  are  greatly  reduced.  A  schematic  for 
the  voltage  amplifier  is  shown  in  Figure  A-4. 

Operational  amplifiers  IC11  and  IC12  serve  as  very  high  impedance  input  buffers  to  avoid 
common-mode-rejection-ratio  (CMRR)  degradation  resulting  from  impedance  mismatch 
between  the  differential  inputs.  Since  the  inputs  to  these  buffers  are  DC  coupled,  the  input 
bias  depends  on  the  signal  source.  Therefore,  high  value  resistors  are  connected  from  each 
input  to  ground  to  provide  DC  return  paths  to  the  amplifier's  common  reference  (Bun’-Brown 
Corporation  1987). 
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Figure  A-4.  Schematic  of  Differential  Voltaae  Amplifier. 
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APPENDIX  B: 

POSITION  DETECTOR  SYSTEM  SOFTWARE 
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B.1  Interpolation  Routine.  The  following  Applesoft  routine  is  designed  to  interpolate 
position  sensor  calibration  data  to  a  full  40  by  40  array  of  points.  Experimental  data  for 
system  calibration  is  not  measured  at  all  1600  positions  of  a  40  by  40  grid  to  avoid  the  tedium 
that  such  a  process  would  involve.  However,  in  order  to  achieve  maximum  resolution  for  the 
corrected  position  value  it  is  necessary  to  know  the  raw  data  values  for  a  40  by  40  grid  of 
points. 


100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 


D$>  CHR$(4) 

DIM  XD(20.20).YD(20.20) 

DIM  XP(20).YP(20).XE(40.40).YE(40.40) 

DIM  XC(20.40).YC(20.40) 

TEXT :  HOME 
VTAB8 

PRINT  THIS  PROGRAM  IS  DESIGNED  TO  READ  IN¬ 
PRINT 

PRINT  THE  RAW  DATA  FROM  AN  EXISTING  FILE- 
PRINT 

PRINT  -AND  INTERPOLATE  THESE  DATA  VALUES  TO  A- 
PRINT 

PRINT  -GRID  OF  DIMENSION  40  X  40- 
PRINT  :  PRINT :  PRINT 

PRINT  -SEE  KOTTKE  IF  YOU  HAVE  ANY  QUESTIONS.- 
VTAB24 

INPUT  -ENTER  A  RETURN  TO  CONTINUE:  -.A$ 

HOME 
VTAB  12 

PRINT  -ENTER  THE  NAME  OF  THE  EXISTING  - 
PRINT  -CALIBRATION  DATA  FILE;- 
PRINT 
INPUT  F1$ 

HOME 
VTAB  12 

PRINT  -ENTER  THE  NAME  OF  THE  INTERPOLATED- 
PRINT  -FILE  TO  BE  CREATED:- 
PRINT 
INPUT  F2$ 

REM  ******^***^*****^**^*^**** 

REM  INPUTTING  CALIBRATION  DATA 

HOME 

PRINT  -INPUTTING  CALIBRATION  DATA- 
PRINT 

PRINT  D$;-OPEN  -;F1$ 

PRINT  D$;-READ  -;F1$ 

INPUT  N 
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480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 


INPUT  M 
FOR  I  -  1  TO  N 
FOR  J  -  1  TO  N 
INPUT  XD(U) 

INPUT  YD(U) 

PRINT  XD(U).YD(U) 

NEXTJ 
PRINT 
NEXT  I 

PRINT  D$;-CLOSE  ";F1$ 

REM  ******^********************* 

REM  CALCULATING  ACTUAL  POSITIONS 
REM  **—**«••*“•***********«• 

HOME 
VTAB  12 

PRINT  “CALCULATING  RAW  DATA  POSITIONS* 

FOR  I  -  1  TO  N 

XP(I)  -  (I  - 1)  *  39  /  (N  -  1) 

YP(I)  -  (I  - 1)  *  39  /  (N  - 1) 

NEXT  I 
HOME 

REM  ***************************** 

REM  INTERPOLATING  CAUBRATION 
REM  DATA  TO  ALL  POSSIBLE 
REM  POSITIONS. 

REM  ***************************** 

PRINT  “INTERPOLATING  CALIBRATION  DATA: 

FOR  I  -  0  TO  39 
FOR  J  -  0  TO  39 

IF  (I  <  20)  THEN  A  »  1:B  -  16:C  -  1 
IF  (I  >  19)  THEN  A  -  16:B  -  1:C  -  - 1 
FOR  K  -  A  TO  B  STEP  C 
FORL-1  TO  16 

TD  -  ABS  (I  -  XP(K))  +  ABS  (I  -  XP(K  +  1))  + 

ABS  (J  -  YP(L))  +  ABS  (J  -  YP(L  +  1)) 

IF  (TD  <  -  2  *  39  /  (N  - 1)  +  0.00001)  THEN 
HX»K:HY»L;  GOTO  860 
NEXTL 
NEXTK 

PRINT  “LOGIC  ERROR  “:  STOP 

D1  -  SQR  ((I  -  XP(HX))  2  +  (J  -  YP(HY))  2) 

D2  -  SQR  ((I  -  XP(HX  +  1))  2  +  (J  -  YP(HY))  ^  2) 

D3  -  SQR  ((I  -  XP(HX  +  1))  ^  2  +  (J  -  YP(HY  +  1))  *  2) 

D4  -  SQR  ((I  -  XP(HX))  *  2  +  (J  -  YP(HY  +  1))  2) 

IF  (D1  -  0)  THEN  D1  -  0.01 
IF  (D2  «  0)  THEN  D2  »  0.01 
IF  (D3  -  0)  THEN  D3  »  0.01 
IF  (D4 . 0)  THEN  D4  -  0.01 
DT  -  1  /  D1  +  1  /  D2  +  1  /  D3  +  1  /  D4 
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950  XE(I,J)  -  (1  /  D1)  •  XD(HX.HY)  /  DT  (1  /  D2)  * 

XD(HX  +  1.HY)  /  DT  +  (1  /  D3)  •  XD{HX  +  1.HY  +  1)  / 

DT  +  (1  /  D4)  *  XD(HXHY  +  1)  /  DT 
960  YE(I.J)  -  (1  /  D1)  *  YD(HXHY)  /  DT  +  (1  /  D2)  * 

YD(HX  1.HY)  /  DT  +  (1  /  D3)  *  YD(HX  +  1.HY  +  1)  / 

DT  +  (1  /  D4)  *  YD(HXHY  +  1)  /  DT 
970  PRINT  "X-";!;"  Y-"*J;*  XE--;INT  (XE(U)  *  100)  / 100; 

•YE-";  INT  (YE(U)  *  100)  /  100;’  HX-’;HX;"  HY=’;HY 
980  NEXTJ 
990  NEXT  I 

1000  PRINT  :  PRINT  :  PRINT 

1010  PRINT  ’DATA  READY  TO  BE  DUMPED  TO  DISK’:  PRINT 
1020  INPUT  ’ENTER  A  RETURN  TO  CONTINUE:  ’;A$ 

1030  PRINT  D$;’OPEN ’;F2$ 

1040  PRINT  D$;’DELETE ’;F2$ 

1050  PRINT  D$;’OPEN  ’;F2$ 

1060  PRINT  D$;’WRITE ’;F2$ 

1070  FORI-0TO39 
1080  FORJ  =  0TO39 
1090  PRINT  XE(I.J) 

1100  PRINT  YE(I,J) 

1110  NEXTJ 
1120  NEXT  I 

1 130  PRINT  D$;’CLOSE  ’;F2$ 


B,2  Generate  Calibration  Look-Uo  Table  Routine.  The  following  Applesoft  routine  is 
designed  to  generate  a  calibration  look-up  taybie  for  use  by  subsequent  data  correction 
routines.  This  calibration  is  performed  by  a  point-by-point  comparison  of  all  possible  data 
pairs  and  the  previously  obtained  raw  calibration  data. 


100  D$»  CHR$(4) 

110  DIMXA(10).YA(10) 

120  DIM  XE(40.40).YE(40,40) 

130  DIM  XC(40.40),YC(40,40) 

140  TEXT  :  HOME 
150  VTAB8 

160  PRINT  THIS  PROGRAM  IS  DESIGNED  TO  GENERATE" 
170  PRINT 

180  PRINT  ’A  POSITION  DETECTOR  LOOK  UP  TABLE  BY’ 
190  PRINT 

200  PRINT  "DETERMINING  THE  CLOSEST  FIT  BETWEEN" 
210  PRINT 

220  PRINT  "ALL  POSSIBLE  DATA  VALUES  AND  THE" 

230  PRINT 
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240  PRINT  'INTERPOLATED  CAUBRATION  DATA" 

250  PRINT  :  PRINT  :  PRINT 

260  PRINT  "SEE  KOTTKE  IF  YOU  HAVE  ANY  QUESTIONS.' 
270  VTAB24 

280  INPUT  'ENTER  A  RETURN  TO  CONTINUE:  ',A$ 

290  HOME 
300  VTAB  12 

310  PRINT  'ENTER  THE  NAME  OF  THE  EXISTING  ' 

320  PRINT 'INTERPOLATED  DATA  FILE.' 

330  PRINT 
340  INPUT  F1$ 

350  HOME 
360  VTAB  12 

370  PRINT  'ENTER  THE  NAME  OF  THE  LOOK  UP  TEXT 
380  PRINT  'FILE  TO  BE  CREATED:' 

390  PRINT 
400  INPUT  F2$ 

420  REM  INPUTTING  CALIBRATION  DATA 
430  *****************^********** 

440  HOME 

450  PRINT  'INPUTTING  INTERPOLATED  DATA' 

460  PRINT 

470  PRINT  D$;'OPEN 'iFlS 
480  PRINT  D$;'READ';F1$ 

490  FOR  1-0  TO 39 
500  FORJ-0TO39 
510  INPUT  XE(I.J) 

520  INPUT  YE(I.J) 

530  PRINT  XE(I.J).YE(I,J) 

540  NEXTJ 
550  PRINT 
560  NEXT  I 

570  PRINT  D$;'CLOSE  ';F1$ 

590  REM  DETERMINE  BEST  FIT  POINT  BY 
600  REM  POINT. 

610  REM 

620  FORI-0TO39 
630  FORJ-0TO39 
640  LV-  10000000 
650  FORK-0TO39 
660  FORL-0TO39 

670  TV  -  ABS  (I  -  XE(K.L))  +  ABS  (J  -  YE(K.L)) 

680  IF  (TV  <  LV)  THEN  LV  -  TV:NA  -  1  ;XA(NA)  -  K:YA(NA)  =  L 
690  IF  (TV  -  LV)  THEN  NA  -  NA  +  1:XA(NA)  -  K:YA(NA)  =  L 
700  NEXTL 
710  NEXTK 

720  XC(I,J)  -  0:YC(l  J)  »  0 
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730  FOR  LA  -  1  TO  NA 
740  XC(I.J)  -  XC(U)  +  XA{LA) 

750  YC(I.J)  -  YC(U)  +  YA(LA) 

760  NEXT  LA 
770  XC(U)  -  XC(U)  /  NA 
780  YC(I.J)  -  YC(U)  /  NA 
790  IF(XC(I,J)-  INT  (XC(I.J))  <  ^)  THEN 
XC(I.J)-.  INT(XC(I.J)) 

800  IF(XC(I.J)-  INT  (XC(I.J))  >  -  .5)  THEN 
XC(I.J)  -  INT  (XC(I.J))  +  1 
810  IF(YC{I.J)-  INT  (YC(U))  <  5)  THEN 
YC(I.J)-  INT(YC{I.J)) 

820  IF(YC(I.J)-  INT  (YC(I.J))  >  »  .5)  THEN 
YC(I.J)  -  INT  (YC(I.J))  -*•  1 

830  PRINT  -X-M;-  Y--J;-  XC--;(C(I.J);-  YC-";YC(U) 
840  NEXTJ 
850  NEXT  I 

860  PRINT  :  PRINT  :  PRINT 

870  PRINT  "READY  TO  DUMP  LOOK  UP  DATA." 

880  PRINT 

890  INPUT  "ENTER  A  RETURN  TO  CONTINUE:  ";A$ 
900  PRINT  D$;"OPEN  ";F2$ 

910  PRINT  D$;"DELETE  ";F2$ 

920  PRINT  D$;"OPEN  ";F2$ 

930  PRINT  D$;"WRITE  •;F2$ 

940  FORI«0TO39 
950  FORJ-0TO39 
960  PRINT  XC(U) 

970  PRINT  YC(U) 

980  NEXTJ 
990  NEXT  I 

1000  PRINT  D$;"CLOSE  ":F2$ 


B.3  Data  Acquisition.  Correction,  and  Display  Routine.  The  following  Applesoft  routine  is 
designed  to  load  calibration  data  from  a  specified  file  into  particular  address  locations,  acquire 
A/D  converter  raw  data,  determine  the  corrected  position,  and  graphically  display  the  result 
The  actual  acquisition  of  data  from  the  A/D  chip  is  performed  by  a  call  to  a  machine  language 
driver  routine  which  is  the  next  listing  in  this  appendix. 

When  this  program  is  first  run,  the  data  from  the  specified  calibration  look-up  data  file  is 
read  in  and  deposited  at  particular  computer  memory  locations.  Specifically,  the  corrected  X 
position,  X^,  for  the  raw  data  pair  (X„Y,}  is  located  on  the  memory  page  which  is  displaced  by 
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X,  pages  from  the  base  page  30h  at  the  Y^h  memory  position  on  that  page  (h  indicates  a 
hexadecimal  number).  The  correspondng  corrected  Y  position.  Y^  is  located  on  the  same 
memory  page  at  the  Y^h-t^Oh  memory  position.  This  somewhat  obscure  placement  of  the 
calibration  data  in  memory  allows  the  indirect  indexed  addressing  mode  of  the  6502 
microprocessor  chip  to  be  utilized  for  convenient  £md  rapid  access  to  the  corrected  data 
(Oe  Jong  1982).  In  this  addressing  mode,  a  number  in  a  convenient  zero-page  location  (lAL) 
is  added  to  the  number  in  the  Y  register  to  determine  the  low-order  byte  of  the  address  of  the 
operand  (ADL).  The  high-order  byte  of  the  operand  address  (AOH)  is  determined  by  the  sum 
of  the  sequential  zero-page  location  (lAL+l)  and  any  carry  from  the  ADL  summation  (C). 
Symbolically,  the  indirect  indexed  addressing  mode  can  be  represented  as: 

ADL  »  [lALl  +  [Y]  (B-1) 

ADH=.[IAL+11  +  C  (B-2) 

operand  >  [ADHrADL] .  (B-3} 

where  square  brackets  denote  the  contents  of  a  memory  location  or  address  and  the  colon 
denotes  concatenation.  So.  to  determine  the  corrected  X  value,  X^,  from  a  pair  of  raw  data 
values  (X.,Y.),  the  contents  of  lAL  are  preset  to  Oh.  Y.  is  transferred  to  the  Y  register,  and  the 
contents  of  IAL-»-1  is  set  to  Xr-^30h.  A  simple  load  accumulator  call  (LDA)  using  the  indirect 
indexed  addressing  mode  then  returns  with  X^  in  the  accumulator.  Corrected  values  for  Y  are 
obtained  in  a  similar  manner  if  the  contents  of  lAL  are  preset  to  80h. 


100  LOMEM:6*4096 
110  D$a  CHR$(4) 

120  DIM  X(20,40),Y(20.40) 

130  TEXT  :  HOME 
140  VTAB  10 

150  PRINT  "THIS  PROGRAM  IS  DESIGNED  TO  ACQUIRE" 
160  PRINT 

170  PRINT  -DATA  FROM  THE  GRID.  LOOK  UP  THE  " 

180  PRINT 

190  PRINT  -CORRECTED  POSITION  VALUE  AND  - 
200  PRINT 

210  PRINT  -GRAPHICALLY  DISPLAY  THE  RESULTS.- 
220  PRINT  :  PRINT  :  PRINT 
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230  PRINT  "SEE  KOTTKE  IF  YOU  HAVE  ANY  QUESTIONS." 
240  VTAB24 

250  INPUT  "ENTER  A  RETURN  TO  CONTINUE:  ";A$ 

260  HOME 
270  VTAB  12 

280  PRINT  "ENTER  THE  NAME  OF  THE  TEXT  FILE" 

290  PRINT 

300  PRINT  "THAT  CONTAINS  THE  LOOK  UP  TABLE" 

310  PRINT 

320  INPUT  "DATA:  ";F1$ 

330  PRINT  D$;"OPEN  ";F1$ 

340  PRINT  D$;"READ  ";F1$ 

350  FOR  1-0  TO  19 
360  PRINT  I 
370  FORJ-0TO39 
380  INPUT  X(U) 

390  INPUT  Y(I,J) 

400  NEXTJ 

410  NEXT  I 

420  FORI-0TO19 

430  PRINT  I 

440  FORJ-0TO39 

450  POKE  3  *  4096  +  I  *  256  -t-  J.  X(I.J) 

460  POKE  3  •  4096  +  I  *  256  +  8  *  16  +  J,Y(I,J) 

470  NEXTJ 

480  NEXT  I 

490  FOR  I  =  20  TO  39 

500  PRINT  I 

510  FORJ  =  0TO39 

520  INPUT  X(l  -  20.J) 

530  INPUT  Y(l  -  20,J) 

540  NEXTJ 

550  NEXT  I 

560  FOR  I  =  20  TO  39 

570  PRINT  I 

580  FORJ  =  0TO39 

590  POKE  3  •  4096  +  I  *  256  +  J,  X(l  -  20.J) 

600  POKE  3  *  4096  +  I  *  256  +  8  *  16  +  J,Y(I  -  20.J) 

610  NEXTJ 
620  NEXT  I 

630  PRINT  D$:"CLOSE  ";F1$ 

640  HOME 
650  VTAB  12 

660  PRINT  "LOADING  MACHINE  CODE  ROUTINE." 

670  PRINT  D$;"BLOAD  DETECTOR  MACH  NOPLOT, A$2000" 
680  TEXT:  HOME 
690  VTAB  12 

700  PRINT"  DETECTOR  IS  ARMED" 

710  REM  POKING  S&H  DELAY  AND  JUMP  VECTOR 
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720  POKE  250,2 
730  POKE  10,76 
740  POKE  11,0 
750  POKE  12,2  *16 
760  A  -  USR(O) 

770  XRAW  .  INT  (  PEEK  (255)  *  39  /  255) 

780  YRAW-  INT  (  PEEK  (254)  *  39  /  255) 

790  XC-  PEEK  (3  *  4096  +  XRAW  *  256  +  YRAW) 

800  YC  -  PEEK  (3  *  4096  +  XRAW  *  256  +  8  *  16  +  YRAW) 
810  HUN  0,39  AT  YC 
820  VLIN0,39ATXC 
830  VTAB22 

840  PRINT -RAW  X  = -pCRAW;"  CORRECTED  X  -  •;XC 
850  PRINT  "RAW  Y  »  ",YRAW;*  CORRECTED  Y  -  ";YCr 
860  FOR  T  -  1  TO  1500:  NEXT  T 
870  HOME 
880  GOTO  680 


B.4  Machine  Language  A/D  Driver.  This  6502  machine  language  routine  is  used  to 
acquire  data  from  the  A/D  converter  chip  when  signalled  by  an  active  low  on  the  nonmaskable 
interrupt  line  (NMI).  The  acquired  data  is  stored  in  zero-page  memory  locations  which  can  be 
accessed  the  controlling  high  level  routine. 


1 

ORG  $2000  ;RUNADDR 

2* 

3 

JSR  $FF4A  ;SAVEREG 

4* 

5 

LDA  #$4C 

LOAD  NMI 

6 

STA  $03FB 

JUMP 

7 

LDA  #<NMI 

COMMAND 

8 

STA  $03FC 

AND 

9 

LDA  #>NMI  •,ADDRESS 

10 

STA  $03FD 

11  * 

12 

LDA  #$FF  ;SETDONE 

13 

STA  $FB  ;FLAG 

14* 

15  LOOP 

LDA  #$FF 

COUNTER 

16 

JSR  $FCA8 

WAIT 

17 

STA  $C030 

SPEAKER 

18 

LDY  $FB 

GET  FLAG 

19 

BNE  LOOP 

DONE  ?? 

20 

JMP  STOP 

TO  EXIT 

21  * 
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22  NMI 

LDA  #<PUNT 

;NON  NMI 

23 

STA  $03FC 

;INTERUPT 

24 

LDA  #>PUNT 

;ROUTINE 

25 

STA  $03FD 

‘ADDRESS 

26* 

27 

LDX  $FA 

;SAMPLE 

28  DELAY 

DEX 

HOLD 

29 

BNE  DELAY 

;DELAY 

30* 

31 

LDA  #$00 

;SET  ADO 

32 

STA  $C0C4 

;OHANNEL 

33* 

34 

LDA  $C0C1 

;START  OON 

35* 

36 

NOP 

;DELAY 

37 

NOP 

;TO  ALLOW 

38 

NOP 

;THE  EOO 

39 

NOP 

40 

NOP 

;D7  UNE 

41 

NOP 

;TO  DROP 

42 

NOP 

43TESTX 

LDA  $C0C2 

;OONVERSN 

44 

BPL  TESTX 

;DONE  ?? 

45* 

46 

LDA  $0003 

;GET  DATA 

47* 

48 

STA  $FF 

;STORE  X 

49* 

50 

LDA  #$01 

;SET  ADO 

51 

STA  $0004 

;OHANNEL 

52* 

53 

LDA  $0001 

;START  OON 

54* 

55 

NOP 

;DELAY 

56 

NOP 

;TO  ALLOW 

57 

NOP 

;THE  EOO 

58 

NOP 

;D7  LINE 

59 

NOP 

60 

NOP 

;TO  DROP 

61 

NOP 

62* 

63  TESTY 

LDA  $0002 

;OONVERSN 

64 

BPL  TESTY 

;DONE  ?? 

65* 

66 

LDA  $0003 

;GET  DATA 

67* 

68 

STA  $FE 

;STORE  Y 

69* 

70 

LDA  #$00 

:DONE 
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71 

STA  $FB 

;FLAQ 

72* 

73 

RTI 

:END  NMI 

74* 

75  STOP 

JSR  $F836 

;OLSORN 

76* 

77 

LDA  #$0F 

;SET 

78 

JSR  $F864 

;OOLOR 

79* 

80 

STA  $0054 

;PRI  PAGE 

81 

STA  $0056 

;LOW  RES 

82 

STA  $0053 

;TXT&GR 

83 

STA  $0050 

;GO  GRAPH 

84* 

85 

JSR  $FF3F 

;REP  REG 

86* 

87 

RTS 

;GO  BASIO 

88  * 

89  PUNT 

RTI 

;NON  NMI 
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